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We study the fidelity decay and its freeze for an initial coherent state of two-mode Bose-Einstein 
condensates in the Fock regime considering a Bose-Hubbard model that includes two-particle tunnel¬ 
ing terms. By using linear-response theory we find scaling properties of the fidelity as a function of 
the particle number that prove the existence of two-particle mode-exchange when a non-degeneracy 
condition is fulfilled. Tuning the energy difference of the two modes serves to distinguish the presence 
of two-particle mode-exchange terms through the appearance of certain singularities. We present 
numerical calculations that illustrate our findings, and propose exploiting a Feshbach resonance to 
verify experimentally our predictions. 
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I. INTRODUCTION 

The Bose-Hubbard model became a workhorse to de¬ 
scribe interactions of ultracold bosonic gases trapped by 
neighboring potentials after its striking success with the 
Mott insulator-superfluid transition mm- Its most sim¬ 
ple physical realization, when only two bosonic states 
can be occupied [3], is experimentally obtained by trap¬ 
ping the condensate in a double-well potential [1]. This 
system is interesting because it is the simplest scheme 
for atom interferometry. In addition to interference phe¬ 
nomena [Si, it also exhibits quantum tunneling and 
self trapping effects m i such as Josephson oscilla¬ 
tions [S]. It has even been used to produce and study 
many-particle entanglement |10j and dynamically gen¬ 
erate spin-squeezed states m- Alternative methods 
to optical lattices have been demonstrated by splitting 
a single-component Bose-Einstein condensate (BEC) on 
atom chips either with pure DC magnetic fields [12] or 
by dressing static fields with RE potentials [TSHT^ as 
proposed in [TBj. Understanding the effects originated 
on inter-atomic collisions has already been exploited to 
overpass the classical limit in atom interferometers [6|, 
for example. Here we employ an extended Bose-Hubbard 
Hamiltonian to increase the possibilities in this direction. 

Several interchange terms arise in the Bose-Hubbard 
model in a double-well potential where only the low¬ 
est level in each well is populated and the correspond¬ 
ing wave functions have a small overlap [H [T^. In 
particular, two terms accounting for two-particle mode- 
exchange processes appear in the derivation of the Hamil¬ 
tonian. These terms are often neglected assuming that 
two-particle processes are rare for diluted ultracold gases. 
Yet, Ref. [T7] points out that there is a better agreement 
with the experimental results when these processes are 
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included. In this paper we probe the relevance of these 
terms by studying dynamical properties linked with two- 
particle tunneling processes. We consider the dynamical 
stability of the quantum evolution under small system 
perturbations for the two-mode Bose-Hubbard model us¬ 
ing the fidelity or Loschmidt echo [T51120j . whose decay 
has been studied for different parameter ranges and types 
of perturbations in the Bose-Hubbard model . 

Prosen and Znidaric noticed that fidelity stops de¬ 
caying, staying essentially constant (modulo some os¬ 
cillations) for relatively long times whenever the time- 
averaged expectation value of the perturbation van¬ 
ishes [24]. This phenomenon is called fidelity freeze. It 
was later shown that symmetries can also induce this be¬ 
havior if the diagonal matrix elements of the perturbation 
vanish, e.g. when the perturbation is not invariant under 
the time-reversal symmetry [25]. Note that the freeze of 
fidelity was actually observed in simulations for bosonic 
and fermionic many-body systems [H] but it was 
attributed to the non-linearities introduced by the inter¬ 
actions between the particles. Our purpose is to draw 
attention to this phenomenon and exploit it, within the 
context of the Bose-Hubbard model. 


The paper is organized as follows. In Sect. [TTJ from 
a generalized two-mode Bose-Hubbard Hamiltonian we 
derive a the fidelity freeze Fpr associated to an initial 
macroscopic trial state m- In Sect, [hi] we show analyt¬ 
ically and numerically that the scaling properties of the 
fidelity freeze display a transition in terms of the number 
of particles if the interaction includes two-particle mode- 
exchange terms. In addition, when certain degeneracy 
condition is fulfilled by tuning the energy difference of 
the two levels, the fidelity freeze tends abruptly to zero. 
This yields insight into many-body tunneling processes 
and provides a method to calibrate the system to en¬ 
hance the fidelity freeze. In Sect. |IV| we summarize our 
results and address the possibility to access the parame¬ 
ter range of interest in this paper considering ®^Rb and 
85Rb. 
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II. FIDELITY FREEZE FOR THE TWO-MODE 
BOSE-HUBBARD MODEL 

The fidelity amplitude is the overlap of the time- 
evolution of an initial state under a reference interaction 
Hq with the evolution of the same initial state under a 
slightly different Hamiltonian H = Hq -|- AH [20] : 

f{t) = {^o\Uo{-t)U^{mo). ( 1 ) 

Here, |'I'o) is the initial state under consideration, 
Uo{t) = Texp[—iHQt/h] is the (time-ordered) unitary 
time-evolution associated to the reference Hamiltonian, 
U\(t) is the corresponding time-evolution of the per¬ 
turbed Hamiltonian, and the perturbation strength is de¬ 
noted formally by A. The modulus squared of the fidelity 
amplitude, F{t) = |/(t)p, is the fidelity or Loschmidt 
echo na IS]. Clearly, F{t) is a measure of the sensi¬ 
tivity of the time evolution of j'f'o) to system perturba¬ 
tions. Another interpretation is that of an echo: j'l'o) 
evolves under Hq up to time t, when the system is sud¬ 
denly reversed with respect to time, and then evolves 
under the action of H ; the Loschmidt echo compares the 
whole evolution with the initial state, thus quantifying 
the degree of irreversibility of the system. The operator 
M\{t) = Ux{—t)UQ{t) is referred as the echo operator. 

We consider the generalized Bose-Hubbard model 
Hbr = Hq + V defined by 

Hq = eihi -I- 62^2 + ^ [ni{ni - 1) + n2{n2 - 1 )], ( 2 ) 
H = - {h\b2 + h\k) - ^ [{blf bj + {blr &?] ■ (3) 

As usual, hi = blbi {i = 1,2) is the particle number op¬ 
erator of the zth-mode, with SJ and hi the corresponding 
bosonic creation and annihilation operators, respectively. 
The single-particle energies of each mode are denoted by 
Ci, [/ is the two-particle on-site interaction, Ji is the en¬ 
ergy of the usual (one-particle) Josephson tunneling or 
mode-exchange term, and J 2 is the energy associated to 
two-particle tunneling processes that we probe here. The 
total number of particles n = ni-\-n 2 is a conserved quan¬ 
tity; fixing n, the Hilbert-space dimension is simply n-l-1. 
The Hamiltonian fLen defined through Eqs. and ([^ 
is a generalization of the usual two-mode approximation 
used to describe the bosonic Josephson junction 
Below we use Ji and J 2 as perturbation parameters, re¬ 
placing A in Eq. 0. 

We are interested in the so-called Fock regime {U >> 
Ji > J 2 ) because the fidelity freeze can be observed there. 
We choose the Fock (occupation-number) basis defined 
by 1 ^ 1 ,^ 2 ) = {bii\b^ 2 \)-^/\b\Y^{b\Y-\Q), where |0) is 
the vacuum state; since n = /ii -I- /J ,2 is conserved, we 
use the short-hand notation |/ii) = |/ii,/i 2 ). By defini¬ 
tion, Hq is diagonal in the Fock basis and V has van¬ 
ishing diagonal matrix elements. Then, considering Hq 
as the reference interaction and V as the perturbation 


or residual interaction, the conditions to observe the fi¬ 
delity freeze are fulfilled [25]. The unperturbed spec¬ 
trum is given by = Eq + /i(ei — e 2 — C/n) -|- with 
Eq = € 2 n + Un{n — I)/2, where /r = 0,... n -|- 1 labels 
the Fock states by mode occupation; notice the parabolic 
shape of in terms of /r for non-vanishing U. As it is 
often done we use the Heisenberg time tn = 2-Kh/d as 
the unit of time, where d is the average level spacing of 
the unperturbed Hamiltonian Hq. 

We compute the fidelity decay by noting that M\{t) is 
the time-evolution propagator associated with the time- 
dependent Hamiltonian Vi{t) = Uo{—t)VUo{t) in the in¬ 
teraction picture [55]. We use Dyson’s series on the per¬ 
turbation parameters Jr (r = 1, 2) truncated to the sec¬ 
ond order |28l [29] . This approach is called linear response 
theory. 

We write the fidelity amplitude as /(t) = 1 -f- /i -I- /2 + 
0{Jf), where the first- and second-order corrections (in 
both Ji and J 2 ) read 

A = E E (4) 

= E E ^2 [t; . (5) 

r,s 

Here, the matrix elements of the perturbation in the in¬ 
teraction picture are exp[iH^^i,t] 

with = Ef^ — Fy, and are the expansion coef¬ 

ficients of the initial state in the Fock basis. In Eqs. 0 
and ([^, greek letters represent the basis states and 
r,s = 1,2 stand for the one- or two-particle tunneling 
terms of Hbk- These matrix elements read 

( 6 ) 

where = [((() To second order in the pertur¬ 

bations, the fidelity is 

F(t) = l + 25R(/i) + 2 5R(/2) + |/i|2. (7) 

The time dependence of Eqs. 0 and (1^ appears in the 
(time-ordered) integrals Ip[t; Hi,..., Hp]7where p stands 
for the order in the Dyson’s series. These integrals can 
be expressed recursively as 

Ip+i [t; Hi,..., Hp+i] = 

— ^ J' dti exp[iHiti] Ip[ti, H 2 ,..., Hp^i], (8) 

where Io[t] = 1 defines the initial value of the recur¬ 
sion. These integrals produce terms that oscillate in 
time as long as the frequencies Hp appearing in the 
exponentials do not vanish, i.e. when the unperturbed 
spectrum is non-degenerate. Yet, certain frequency com¬ 
binations may vanish and yield secular terms which grow 
at least linearly in time. We assume that the unper¬ 
turbed spectrum is non-degenerate, which can be assured 
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by choosing properly the energy difference of the two 
modes Ae = €2 — ei- Then, without the secular contribu¬ 
tions, to second-order the fidelity displays quasi-periodic 
oscillations in time; this is the freeze of the fidelity. The 
freeze of the fidelity lasts as long as the second-order ap¬ 
proximation is valid; eventually, higher-order contribu¬ 
tions dominate the evolution and secular terms appear 
that destroy the freeze of the fidelity. 

Equation Q is valid for any initial state. We consider 
as the initial state a normalized macroscopic trial state 
of the form m 

i^o} = (abl+/^e^^blriO} 

/ \ 1/2 

= Zl(j (9) 


that Xi[t;r2] = (1 — exp(iflt))/hfi, and X 2 [^; fli, 1 ^ 2 ] = 
+ 0 , 2 ])/Kl 2 ^ for and ^2 non-zero. 

These conditions are fulfilled by the assumption of a non¬ 
degenerate spectrum. Yet, for fli - 1-^2 = 0 a secular term 
is obtained for /2 which has the form —it/h. This term 
does not affect the fidelity according to Eq. Q , since it is 
purely imaginary. Then, to second-order in the tunneling 
rates, the time-dependence of the fidelity is at most quasi- 
periodic, hence fidelity exhibits a freeze. The time during 
which the fidelity freeze lasts scales as the inverse of the 
perturbation and the inverse of n; see P5] . 

In order to obtain the fidelity freeze we extract the 
time-independent contributions of the integrals Q in the 
expression for 3?(/i), 3?(/2) and |/ip. Including the de¬ 
pendency of on the phase (j), cf. Eq. ([^, we obtain 


This initial state is coherent HZ]; with a = 
and j3 = (n 2 /n)^/^, it corresponds to the mean-field state 
having ni particles in the first mode and n 2 = n — ni in 
the second one. 

Inserting Eqs. ^ and (|^ in Q and ([^, we obtain 








( 10 ) 


/2 =J 2 jrJs 

r,s 

X (I (I 1 

^ y^,n-p'-’p,p.-r\yp,n-v'Jv,p-s ' Un-p,u^^',p-s) 

~^9ri-p,p^P,P+r{9p,n-v^i',p-s + 9n-p,i^'^i',P+s) ■ (H) 

The time dependence can be further described by noting 


3? 


3? 


3? 




cos {p(j)) 


hn 


e^P^l 2 [t-M,,n 2 ] 




cos(pcj)) / 1 
^2^2 vfii; 
C0s(p(/>) 


1 — ^ 1,-2 
rii -f f22 

(1 -I- 5i,2)- 


( 12 ) 

(13) 

(14) 


h‘^nin 2 

Here, the right-hand side of these expressions are the 
time-independent contributions, where p is an integer 
related to the indexes of the Fock states involved, and 
we have used the Kronecker-delta 81-2 to indicate that 
the frequencies satisfy Hi = —H 2 (indexes are reversed), 
and Si _2 to denote that Hi = H 2 (indexes are the iden¬ 
tical). Note that in Eq. (13) the secular term related to 


Hi -I- H 2 = 0 is not included due to the (5i._2. 

Inserting Eqs. (10) and (11) into (^ and using the 


time-independent contributions, Eqs. ([r2|) to (|14|), we ob¬ 
tain 


Fpr = 1 -k 2 ^ J, 


cos{r(j)) 


hfi 




"h 2 ^ ^ JrJs 


cos((Tr T s)^) 


1 - 6 . 




/j.±r,/.A±r±s 


.H 


p,p±r 


H 




p,p±r±s 


^ ^ \Al^A^±r\\AlA^±s\Q\^'''^Q^^"r(^os[{Tr±s)(j)) 

^-(l+<)±r.±«). 




(15) 


The coefficients and = gnlp,p+r 

are introduced to have a more compact expression. The 
signs of r and s are independent and correspond to the 
distinct possibilities imposed by the Kronecker deltas 
that appear in Eqs. (10) and Equation ( [I^ is a 

central result of this paper. 


In Fig. we show an example of the decay of fidelity 
for a coherent state with ni = 712 = n/2 and </) = 7 r /4 
obtained numerically. The figure illustrates the oscilla¬ 
tions during the freeze of the fidelity, the eventual decay. 


and the value obtained from Eq. ( |T^ for the freeze of 
the fidelity (horizontal green line). Time is measured in 
Heisenberg-time units tn- The parameters of the model 
are 17 = 1, Ji = 10"®, J 2 = 10"®, ei = 0.76, e 2 = 0.93 
and n = 128; the values of assure the non-degeneracy 
of the spectrum of Hq (see below). These parameters 
have been chosen to simplify the numerics; other values 
display qualitative similar behavior as long as we are in 
the Fock regime. In the inset we display the result consid¬ 
ering the second-order expansion ([^; the value of Fp^ is 
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III. SCALING PROPERTIES OF THE 
FIDELITY FREEZE 


We address now the scaling of i^pr in terms of the num¬ 
ber of particles. An estimate of the scaling properties is 
obtained considering the maximum contribution of the 
n-dependent terms in Eq. (15). This follows from a Fock 
state that we write as /i = An, and then use Stirling’s for¬ 
mula for large n. It can be shown that \A*^A^±r\ 

and ~ n® for A = . The scaling laws of the time- 

independent contributions thus read 5i(/i) ~ 

3?(/2) ~ J> 2 r-i /2 3fJ(|/i|2) ^ J2n2’'-1. Hence, the 

asymptotic dominating contribution for the fidelity freeze 
scales as 


I-Fr. ~ (jg) 


FIG. 1: (Color online) Time-dependence (in Heisenberg time 
units) of the fidelity obtained numerically for the extended 
Bose-Hubbard model in log-log scale, using an initial coher¬ 
ent state Eq. Q with ni = n 2 = n/2 = 64 and (j) = 7r/4. 
The parameters for the Hamiltonian are U = 1, Ji = 10“®, 
J 2 = 10“®, £1 = 0.76 and £2 = 0.93. The horizontal line cor¬ 
responds to the value of the fidelity freeze Fpr obtained from 
Eq. (151. In the inset we present the result of the second-order 
linear-response theory. 



logid") 

FIG. 2: (Color online) Log-log plot of 1 — Fpr as a function 
of the number of particles, for an initial symmetric coherent 
state Til = n 2 = n/2. Triangles (green, dotted line) cor¬ 
respond to the parameters Ji = 10 “® and J 2 = 0 , squares 
(blue, dashed line) to Ji = 0 and J 2 = 10 “®, and circles (red, 
continuous line) to Ji = 10 “® and J 2 = 10 “®; the remaining 
parameters are those of Fig. 


an average over the quasi-periodic oscillations that take 
place during the freeze. 


This result predicts a different scaling for each of the tun¬ 
neling terms Ji and J2. Thus, Fpr exhibits a transition 
from a behavior dominated by Ji to a regime where J 2 
dominates, around n ^ J1/J2. Figure]^ is the numeri¬ 
cal confirmation of this statement. The data points were 
obtained numerically from the time series (cf. Fig. [^, 
using the maxima of the quasi-periodic oscillations of 
1 — F{t) during the freeze; these values underestimate 
the theoretical expectation for Fpr. Fitting the data to 
straight-lines when either J 2 or Ji are absent yields the 
slopes 1.52 and 3.54, respectively. These values are in 
excellent agreement with the 3/2 and 7/2 predicted by 
Eq. (16), thus showing that the scaling properties of the 
fidelity freeze in terms of n probe the presence of two- 
particle tunneling processes. Equation (16) remains valid 
for small departures from the symmetric initial coherent 
states, i.e., A ~ 

An import ant assumpt on that we made in the deriva¬ 
tion of Eqs. (15) and (16) is that the spectrum of Hq is 
non-degenerate, which can be fulfilled by tuning Ae, the 
energy difference of the two modes. As we approach a de¬ 
generacy, the appearance of secular terms makes Eq. (15) 
no longer valid. This can be exploited to probe the rele¬ 
vance of two-mode exchange processes. 

To clarify this idea we consider the Fock state /iq = 
[(n + Ae/C/)/2] whose energy is the minimum of the 
spectrum of i7o, where [cc] is the round-to-nearest inte¬ 
ger function; note that this minimum corresponds to the 
parabolic shape of induced by a non-vanishing U. As¬ 

suming that n is even for concreteness, it can be shown 
that AejU = 0 implies that F^o-i = A^o+Ij meaning 
that the Fock states /ig — 1 and /ig + 1 are degener¬ 
ate; these states are coupled by a two-mode tunneling 
term. The same holds for AejU = 2, though the ac¬ 
tual value of /Tg has changed. For Ae/U = 1 we have 
= F^(,_|_i, i.e. the ground state is degenerate, which 
also holds for AejU = 3; in this case, the states are cou¬ 
pled by a one-particle tunneling term. Then, by tuning 
the single-particle energies, as we approach AejU = ±1 
or ±3, a peak in log]^g(l — F-pA) develops indicating that 
the perturbation does contain a one-particle tunneling 
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FIG. 3: (Color online) Behavior of log]^Q(l — Fpr) as a func¬ 
tion of the energy difference between modes Ae scaled by 
the two-particle interaction coefficient U of the Bose-Hubbard 
model, Eq. The 3D plot depict the appearance of a peak 
around AejU — 0,2 which becomes noticeable as the particle 
number n increases. 

term; likewise, a peak at AejU = 0 and ±2 appears if 
there are two-particle tunneling processes. This is illus¬ 
trated in Fig. which depicts logiQ(l — Fpr) in terms of 
AejU for various even values of n. Note that the narrow 
peaks at AejU = 0, ±2, the signature of the two-particle 
tunneling, grow for increasing values of n. For odd val¬ 
ues of n the same argument applies, exchanging only the 
location of the peaks. Thus, by increasing n, the peaks 
associated with the two-particle tunneling processes be¬ 
come comparable to those associated to the one-particle 
tunneling processes; for big enough n the distance be¬ 
tween prominent neighboring peaks is halved. This result 
means that the fidelity freeze Fpr can also be maximized 
by tuning AejU. 


IV. SUMMARY AND OUTLOOK 

Summarizing, we have found that the fidelity freeze 
from an initial symmetric coherent state is a sensitive 
quantity to two particle mode-exchange processes in the 
Bose-Hubbard model. This sensitivity can be controlled 
with two experimental parameters: the total atom num¬ 
ber n and the energy difference between modes Ae/U . In 
terms of n, the fidelity freeze displays a transition from a 
regime dominated by the one-particle exchange term, for 
small particle numbers, to the dominance of two-particle 
tunneling processes when n is large enough (n ~ J\l 
There, the fidelity freeze can also be maximized by tuning 
Ae/U. 

Our findings hold in the Fock regime of a double well 
potential, i.e. for J^/U <C 1. Their test would face two 
technological challenges: measuring fidelity and produc¬ 



FIG. 4: Gomputed tunneling parameters of the two-mode 
BEG model as a function of the potential-barrier height, 
where the trap considered corresponds to that of Ref. [8]. 
The red curves correspond to ®^Rb, and the blue, green and 
purple to ®®Rb with s-wave scattering lengths of 20a, 50a and 
100a, respectively, where a — 5.32nm is the s-wave scattering 
length of ®^Rb. The upper curves correspond to Ji/U and 
the lower to J 2 /U. 

ing a BEG with an adequate atom number and confin¬ 
ing geometry. Measuring fidelity is not a simple task 
though it has been achieved in NMR polarization echo- 
spin experiments |.30L I31j and in periodically-kicked cold 
atoms [32] • Echo spectroscopy experiments in cold atoms 
by Andersen et al. [33] demonstrated measurements of a 
quantum fidelity defined differently; theoretical aspects 
of that definition are discussed in Ref. [31]. Fidelity has 
not yet been measured in two-mode Bose-Einstein con¬ 
densates, though it has been proposed in the context of 
cold optical lattices Ref. [85] . 

Regarding what is an adequate BEG species to probe 
our results, special care is required for having a long 
enough confinement with J 2 big enough, so the effects 
addressed here can be observed experimentally. As an 
example, consider ®^Rb using the same trap frequencies 
and separation of the wells as in Ref. [5], and vary the 
potential height in order to reach the Fock regime. In 
this case, the Heisenberg time is ^ 5 — 7 ms for 1500 
bosons. The duration of the fidelity freeze is x 10^ 
(Fig. [I]), which is well inside normal experimental times 
for Bose-Einstein condensates under ultra-high vacuum. 

In Fig. 1 ^ we present the parameters Ji/U (upper 
curves) and J2/U (lower curve) using the formulae of 
Ref. [17], as a function of the barrier height Vq for n = 
1500 bosons. These results were obtained by integrating 
the non-polynomial non-linear Schrodinger equation |36j 
using the standard split-slit Fourier method [371[3E]. The 
results show that for ®^Rb (red curves), the two-particle 
tunneling parameter J 2 /U is perhaps too small to yield 
any measurable signal. For instance, for Vq ~ I500Hz, 
where Ji/U ~ 10“^ — 10“^, we obtain J2/U ^ 10“^^. 
Increasing further the potential height makes J2/U to 
decrease even further. 

A more promising possibility is to consider other mech- 
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anisms that increase J 2 /U, e.g. approaching a Feshbach 
resonance; an obvious candidate is ®®Rb [39]. Consid¬ 
ering the same parameters for the trap used above, the 
blue, green and purple curves in Fig. [^ correspond to 
s-wave scattering length values of 20a, 50a and 100a, re¬ 
spectively. Here we use a = 5.32 nm, the s-save scattering 
length of ®^Rb, as unit to ease the comparisons. The re¬ 
sults show that larger values of J 2 /U are obtained and, 
in that sense, may be accessible to experimental observa¬ 
tion. Yet, we note that increasing the s-wave scattering 
length by such amounts leads to three-body collisions 
which have not been taken into account in our calcula¬ 
tions. 

Double-well potentials are now a common scenario for 
atom interferometry with matter waves. Our findings 


could be useful to study coherence and decoherence ef¬ 
fects in this context. For instance, they could help 
to minimize decoherence on interferometers using dense 
atomic clouds, where non-linearities due to collisions can 
be exploited to improve their accuracy |S] , or give rise to 
optimal methods for analyzing the interference fringes 
imprinted by small energy differences between matter 
waves m- 
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